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We show that polarization entangled photons at x-ray energies can be generated via spontaneous 
parametric down conversion. Each of the four Bell states can be generated by choosing the angle of 
incidence and polarization of the pumping beam. 
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Spontaneous parametric down-conversion at hard x- 
ray energies was first proposed by Freund and Levinc in 
1969 [l[, and first demonstrated experimentally by Eisen- 
bcrgcr and McCall about two years latter [2|. In that 
work a hard x-ray tube was used as the pump, and coin- 
cidence counts at the signal and idler were measured at 
the rate of a few counts per hour. The first experiment 
using a synchrotron was done in 1997 by Yoda ct al. Q 
with a counting rate of 6 counts per hour. During the 
years 1998-2003, Adams and collaborators conducted a 
series of experiments and improved the coincidence count 
rate to about 1 count per 13 seconds Ql- Recent and 
expected improvements in brilliance and beam quality 
of synchrotron x-ray sources, together with new facilities 
such as the x-ray free electron laser and energy recovering 
linacs offer the possibility of extending the concepts 
of quantum optics as developed in the visible portion of 
the electromagnetic spectrum [H to x-ray wavelengths. 

As a step toward this extension, this Letter describes 
a method for generating polarization entangled photons, 
in pure Bell states, at x-ray wavelengths. The technique 
is straight-forward and makes use of the selection rules 
that are associated with a phase matched plasma-like x- 
ray nonlinearity Q ■ In the following paragraphs we will 
show that by choosing the polarization and angle of in- 
cidence of the pumping beam, and working off of the 
degenerate frequency, that each of the four Bell states 
may be generated. A consequence of this work is that, in 
each of the previous x-ray down conversion experiments 
mentioned above, the generated photon pairs were po- 
larization entangled, but in no case were they in a pure 
Bell's state. 

Before proceeding we note two previous suggestions for 
generating entangled photons at x-ray energies. The first 
is a proposal by Schiitzhold et al. who have suggested the 
use of ultra-relativistic electrons accelerated by a strong 
periodic electromagnetic field (for example, a laser or un- 
dulator) to create entangled photon pairs in the multi- 
keV regime @. The second is a proposal by Palffy et al. 
who suggest generating single-photon entangles states by 
control of nuclear forward scattering [3] . 

We start by discussing the nonlinearity. The central 
concept of the nonlinearity at x-ray energies is that, since 
x-ray photons have energies that are large as compared to 
the electron binding energy of low-Z atoms, that the x-ray 



nonlinearity of an element such as diamond may be cal- 
culated by treating all of the electrons in the atom as free 
particles, and therefore treating the nonlinear medium as 
a very dense cold plasma [J, Q . This x-ray nonlinearity is 
of second order so that two frequencies may add or sub- 
tract to generate a third frequency. Three processes that 
are at first glance seemingly different contribute to the 
x-ray nonlinearity. These are: 1) A Lorentz force term 
where the electron velocity caused by an incident electric 
field at frequency wi mixes with the magnetic field of fre- 
quency UJ2 to generate a force and current at frequency 
UJ3 = uji + UJ2- 2) A term that depends on the spatial 
variation of charge density, and 3) A term that depends 
on the spatial variation of velocity. Each of these pro- 
cesses produces a driving current with a k- vector that is 
the sum of the k vectors of the applied fields and of the 
lattice. 

In order to satisfy permutation symmetry and to con- 
serve photons in a three-frequency nonlinear optical pro- 
cess, it is essential that the three processes of the previous 
paragraph all be retained in the calculation of the non- 
linearity Q. It is easily shown that each of the above 
processes, for example the Lorentz force process, does 
not in its own right satisfy permutation symmetry. 

A typical phase matching diagram for parametric down 
conversion is shown in Fig. [TJ With the k- vectors kg, ki, 
kp denoting the k- vectors of the signal, idler, pump fields 
and G denoting the reciprocal lattice vector, the phase 
matching condition for parametric down conversion is 
ks + ki = kp + G. With the unit vectors ej denoting the 
polarization of the respective electric fields, these fields 
are written as Ej{r, t) ~ -f- exp[— — kj ■ r)\ej + c.c. 
Working in the cold plasma approximation, we perturba- 
tively calculate the nonlinear current density Jg (r , t) = 
Ps{r,t)vs{r,t) at the signal frequency 0,0, [lOl to obtain 
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Here q and m are the electron charge and mass, and 
poir) is the electron density in the absence of the pump- 
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FIG. 1: Phase matching diagram for x-ray spontaneous down 
conversion . ks, ki, and kp are the wave vectors of the signal, 
idler, and pump fields. G is the reciprocal lattice vector of 
the diamond crystal. 



ing beam. We substitute the expressions for the electric 
fields into Eq. (jlj and project the nonlinear current den- 
sity against the direction of the signal electric field. We 
assume phase matching [l|, with the reciprocal lattice 
vector G so that the unperturbed electron charge den- 
sity is taken as po{r) = pgexp[iG ■ r\. The envelope of 
the nonlinear current density is then 
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As shown in Fig. [T] we define the scattering plane as the 
plane containing the k- vector of the pumping beam and 
the lattice k-vector, and assume that the k-vectors of 
the signal and idler beams are also in this plane. From 
Eq. ^ we find the following selection rules Q: (1) If 
the polarization of the pump is in the scattering plane, 
the polarizations of the signal and the idler photons must 
both be either in the scattering plane or must both be 
normal to the scattering plane. (2) If the polarization 
of the pump is normal to the scattering plane, then ei- 
ther the signal polarization is in the scattering plane and 
the idler polarization is normal to the scattering plane or 
vise-versa. Polarization entanglement requires that the 
polarization of the idler is uniquely determined by the 
polarization of the signal. For the pump polarized either 
in, or orthogonal to, the scattering plane the down con- 
verted signal and idler photons are therefore entangled. 

The polarization of the entangled photon pairs is de- 
termined by the driving current as described by Eq. (2) 
and is not influenced by the temporal or angular disper- 
sion of the system. To calculate Glauber correlations and 
the biphoton generation rate, we work in the Heisenberg 
picture and write a pair coupled equations for each of the 
biphoton pairs For example, for the entangled state 
|(/) >= {\Hs, Vi > +\Vs,Hi >)/-\/2, we write a pair of cou- 
pled equations for the state \Hs,Vi > and a second set 



of coupled equations for the state \Vs,Hi >. It is critical 
that the dependence on frequency and angle of emission 
of the k-vector mismatch function is the same in each pair 
of coupled equations. Each of the biphoton wave pack- 
ets \Hs,Vi > , or \Vs,Hi > has significant dispersion and 
the temporal and spatial correlations between these pack- 
ets will vary with position and angle. But because the 
k-vector mismatch is independent of polarization and is 
the same for each packet, the polarization correlations are 
determined by the selection rules associated with Eq. (2), 
and do not vary with propagation. 

We next describe how to generate the four maximally 
entangled 2-qubit Bell states. We denote \H > as the po- 
larization of the x-ray electric field in the scattering plane 
(i.e., the plane containing the incident k-vector and the 
lattice k-vector G), and > as the polarization orthogo- 
nal to the scattering plane. With the pump polarization 
in the scattering plane, the polarization of the emitted 
photon pair is 

1^ >^^[Ai6p)\Hs,H, > +Biep)\Vs,V,>] (3) 

Here 9p, is the angle of the pump k vector with regard to 
the atomic planes that are normal to G. The coefficient 
A{9p) is the nonlinear current density when the polariza- 
tion of both signal and idler are in the scattering plane, 
and the coefficient B{9p) is the current density when the 
polarization of both signal and idler are normal to the 
scattering plane. Similarly, when the pump polarization 
is normal to the scattering plane 
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The coefficient C{9p) is the current density when the 
signal is polarized in the scattering plane and the idler 
is polarized normal to the scattering plane. The coef- 
ficient D{9p) is the current density when the signal is 
polarized normal to the scattering plane and the idler 
is polarized in the scattering plane. The quantities 
A{9p), B{9p),C{9p), and D{9p) are real functions of 9p. 
This is different than conventional spontaneous para- 
metric down conversion in the optical regime where the 
equivalent coefficients are generally complex [l2. I3|. 

The implication of Eq. ^ and (g]) is that the proba- 
bilities for generating the states \Hs,Hi > and > 
are |^(0p)p and \B{9p)\'^ respectively; and the probabil- 
ities for generating the states \Hs,Vi > and \Vs,Hi > 
are \G{9p)\'^ and \D{9p)\'^. Consequently, unless either 
\A{9p)\^ = or \B{9p)\^ = 0, Eq. © describes an en- 
tangled state. A maximally entangled state is obtained 
when \A{9p)\ = \B{9p)\. Similarly, Eq. (g]) describes 
an entangled state which is maximally entangled when 
|C(6'p)| — |D(6'p)|. As we will show below, to produce all 
four of the Bell states, it is required that a;,, is not equal 

to UJi. 
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FIG. 2: Polarization entanglement at degeneracy : (a) Square 
of the current density for the states \Hs,Hi > (solid green) 
and \Vs,Vi > ( dashed red), i.e. |^(6'p)p and \B{ep)\'^ as func- 
tion of the pump angle 9p. (b) Square of the current density 
for the states \Hs,Vi > (solid green) and \Vs,Hi > (dashed 
red), i.e. |C(6lp)|^ and \D{dp)\^. The inset in (a) shows the 
intersection of the states \Hs,Hi > and |Vs, Vi >. The inset 
in (b) shows the region containing the two intersection points 
on the left hand side of curves of the states \HB,Vi > and 
\Va,Hi >. A similar region exists on the right hand side. 



The procedure that we use to determine the Bell states 
is to first plot the square of the current density (Fig. 2) 
for each component of the Bell states. The intersections 
of the component curves determine the pump angles at 
which the magnitudes of the components of each Bell 
state are equal. We then, numerically, by simultaneous 
solution of the phase matching and current density equa- 
tions, determine the sign of the components at the in- 
tersection. Consider a specific example: We choose dia- 
mond for the nonlinear medium and use the (111) lattice 
k-vector for phase matching. We take the pump pho- 
ton energy as 25 keV, and first consider the frequency- 
degenerate case where both the signal and idler energies 
are 12.5 keV. We solve the phase matching equations for 
the angles of the signal 9s and idler 9i with regard to the 
atomic planes, and substitute the related electric fields 
into Eq. (2). With the polarization of the pump, signal, 
and idler chosen, the current density is a function of Op 
only. We calculate | Js P for each of the polarization states 



TABLE I: Angles of the k vectors of the pump, signal, and 
idler for producing maximally entangled states at the degen- 
erate frequency. 



Bell's state 


Op 


Os 


Oi 


{\Hs,V^ > +\Vs,Hi >)/V2 


0.1208 


-0.1208 


-0.1208 




0.2434 


-0.2434 


0.2434 




2.89819 


2.89819 


3.385 




3.02079 


3.26239 


3.26239 


i\H,,V, > -\Vs,H^ >)/^/2 


7r 
2 


2.2798 


0.8617 


TABLE II: Angles of the k vectors of the pump, si{ 


;nal, and 


the idler for producing maximally entangled states at 20% off 


of the degenerate frequency. 








Bell's state 


Op 


Os 


0^ 


i\H,,Hi > +\Vs,Vi >)/V2 


1.3252 


0.708018 


2.13134 


i\H,,H, > -\V,,V, >)/V2 


1.15798 


0.512083 


1.8803 


i\Hs,V, > +\Hs,V^ >)/V2 


0.525749 


-0.0629786 


0.84284 


{\Hs,V, > -\H„V, >)/V2 


1.31096 


0.690752 


2.11046 



and plot the results in Fig. [2l Figures 2(a) and 2(b) 
show |A(6'p)|2 and |B(0p)P, and |C(6'p)|2 and \D[ep)\^, 
respectively, all as function of Op. From Fig. [2] (a) we 
see that |A(6'p)p = \B{9p)\'^ only when the nonlinearity 
is zero. Therefore at the degenerate frequency, and with 
the pump polarized in the scattering plane, the gener- 
ated photon pairs are polarization entangled, but a pure 
Bell's state cannot be produced. 

On the other hand, the solutions for |C(6'p)p = 
|D(6'p)p are obtained at finite nonlinearity therefore al- 
lowing the Bell states -i=(|T4,iJj > ±\Hs,Vi >) to 
be generated. As shown in Fig. H^b), the equation 
l^(^p)P = \D{Op)\^ has five solutions; i.e. there are five 
intersections between of the dashed and solid curves. One 
solution is obtained at 6'p = ^ and corresponds to the 
state -^{\Vs,Hi > — \Hs,Vi >). The other four solutions 

correspond to the state ^(|14,-ffi > +\Hs,Vi >). The 
Bell states and the corresponding angles of the pump, 
signal and idler are determined by numerically solving 
the phase matching and current density equations. The 
results are summarized in Table HI 

Next, we analyze a case where the signal frequency is 
20% off degeneracy as is shown in Fig. [31 In this case 
it is possible generate each of the four Bell states. The 
corresponding angles of the pump, signal and idler are 
summarized in Table HIl We note that in contrast to the 
degenerate configuration, when off of degeneracy there is 
only one solution for each of the Bell states. That is, the 
angles of pump, signal and idler fields with regard to the 
atomic planes are uniquely defined by the Bell state. 

To estimate the efficiency for the generation of Bell 
polarization states we solve the coupled Heisenberg- 
Langevin equations for each of the biphoton pairs, and 
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FIG. 3: Polarization entanglement at 20% off of the degen- 
erate frequency (a) Square of the current density for the 
states \Hs,Hi > (solid green) and |K,Vi > ( dashed red), 
i.e. |^(Sp)|^ and \B{9p)\^ as function of 9p. (b) Square of 
the current density for the states \Hs, Vi > (solid green) and 
\Vs,Hi > (dashed red), i.e. |C(0p)p and \D{ep)f as function 
of Op . The inset in (a) shows the intersection of the states 
\H„H, > and \Vs,V, >. 

numerically calculate the generation and coincidence 
count rates . For diamond, for the state {\Hs,Vi > 
— \Hs,Vi >)/\/2 with the signal frequency 20% above the 
degenerate frequency, a crystal length of 2 mm, a pump 
flux of 10^'^ photons/sec (available at the brightest syn- 
chrotron facilities), and detector apertures of 5 mrad x 
5 mrad, the estimated coincidence count rate is about 1 
count per 15 seconds. This count rate is comparable to 
the count rate of previous experiments 0] ■ 

At x-ray wavelengths a polarizer may be constructed 
by Bragg scattering with a Bragg angle of 0^ = 45°. 
When this is the case, the coefRcient for scattering at 29b 
is (ideally) zero when polarized in the plane of incidence 



and unity when polarized perpendicular to the plane of 
incidence. Bragg polarizers with an energy bandwidth of 
AE/E ~ 10~^ might be constructed using mosaic crys- 
tals such as pyrolytic graphite fl^ . 

In summary, this work has described a technique for us- 
ing parametric down conversion at x-ray wavelengths to 
generate each of the four Bell polarization states. When 
off-degenerate this is done by choosing the angle of inci- 
dence of the pumping laser and polarizing it either in, or 
out, of the plane of incidence. When at degeneracy, the 
pump must be polarized out of the plane of incidence, 
and only two of the four Bell states may be obtained. 

The authors thank Jerry Hastings for suggesting the 
x-ray polarizer, as above. This work was supported by 
the U.S. Air Force Office of Scientific Research and the 
U.S. Army Research Office. 
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